Magnetic response of mesoscopic superconducting rings with two order parameters 
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The magnetic response and fluxoid transitions of superconducting aluminum rings of various 
sizes, deposited under conditions likely to generate a layered structure, show good agreement with 
a two-order-parameter Ginzburg-Landau model. For intermediate couplings, we find metastable 
states that have different phase winding numbers around the ring in each of the two order parame- 
ters. Those states, previously theoretically predicted, are analogous to fractional vortices in singly 
connected samples with two-order-parameter superconductivity. Larger coupling locks the relative 
phase so that the two order parameters are only manifest in the temperature dependence of the 
response. With increasing proximitization, this signature gradually disappears. 
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Since the discovery of two-gap superconductivity in 
MgB2 the coexistence of two not-too-strongly cou- 
pled superconducting order parameters (OPs) has mo- 
tivated significant theoretical work. A striking predic- 
tion is the existence of vortices carrying unquantized flux 
0, [H . For non- negligible Josephson coupling, those ex- 
hibit a soliton-shaped phase difference between the two 
OPs 0, HI- Theoretically, such solitons may also form 
when current flow along a wire causes the relative phase 
to unlock, leading to a higher current than in the phase 
locked state Most of this theoretical work is based 
on a two-OP Ginzburg-Landau (GL) theory [f|. In the 
realm of superconductivity, this model applies for both 
Josephson coupled bilayer systems, which are described 
microscopically by a tunneling BCS Hamiltonian, and in- 
trinsic two-gap systems Q ■ A particularly interesting ex- 
ample of the latter is Sr2Ru04. Its p x + ip y OP would im- 
ply zero energy core excitations with non-abelian braid- 
ing statistics, which have been envisioned as a basis for 
topologically protected quantum computation Hi]. A 
similar two-OP model has also been used to describe a 
liquid metallic state of hydrogen [l(| , where the electrons 
and protons would form two independent supcrfluids. 

Compared to the large number of theory papers, there 
is little experimental work on the mesoscopic structure 
of such two-OP systems. In this paper, we report the ob- 
servation of soliton states and other phenomena arising 
from the interplay between two OPs in quasi-lD, super- 
conducting rings consisting of two parallel, Joscphson- 
couplcd aluminum layers. Positioning a scanning SQUID 
microscope 11| over each ring individually enabled mea- 
surements of the current, /, circulating the ring as a func- 
tion of applied flux, $ a , and temperature, T. The ensem- 
ble of magnetic responses includes distinct features that 
cannot be explained by one-OP GL, but can be described 
by numerical solutions of two-OP GL. The inferred cou- 
pling between the two OPs depends on the ring's annu- 
lus width, w, allowing us to study the crossover between 
intermediate and strong coupling regimes. For interme- 



diate coupling, we find metastable states with different 
phase winding numbers for each OP. Those are the ID 
analogue of unquantized vortices 0] , and imply a soliton- 
shaped phase difference 0]. In this regime, multiple tran- 
sition pathways between states give a rich structure of 
hysteretic <& a -I curves. At stronger coupling, the sys- 
tem approaches the Cooper limit of complete proximiti- 
zation [12] : the $ a -7 curves at any temperature can be 
described by one-OP GL, but the existence of two OPs 
is manifest in the temperature dependence of the fitted 
penetration depth, A, and the GL-coherence length, £gl • 

During a single, two-month-long cooldown, we charac- 
terized the magnetic response of 40 different rings with 
eight annulus widths 45 nm < w < 370 nm and radii 
R of 0.5, 0.8, 1.2 and 2 /mi. The rings were fabricated 
on oxidized silicon, using liftoff lithography with PMMA 
resist. The 40 nm thick Al film was deposited by e-beam 
evaporation at a rate of about 1 A/s and a pressure of ap- 
proximately 10~ 6 mBar. During the deposition, the rate 
temporarily dropped to a negligible level for about 10 min 
and subsequently recovered, which most likely caused the 
formation of two superconducting layers separated by an 
AlO x tunneling barrier. A disk with a radius of 2 /itm 
had a T c of 1.6 K, representative of the bulk film. Using 
£o = 1-6 /urn for pure bulk Al [3] and £gl(0) ~ 70 nm 
for our rings, as derived below, we infer a mean free path 
of l e = 1.4£gl(0) 2 /£o = 4 nm. The measured critical 
temperatures T c of the rings ranged from 1.5 to 1.9 K, 
depending only on if. Both the short l e and the large 
T c compared to clean bulk Al indicate small grains and a 
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(likely related) strong effect of oxygen impurities 

Our SQUID sensor has two counterwound pickup loops 
and field coils that are used to apply a local magnetic field 
[l6| . We position one pickup loop over a ring, record time 
traces of the SQUID response while sinusoidally varying 
<& Q at a few Hz, and average hundreds to thousands of 
field sweeps. A background, measured by retracting the 
SQUID from the sample, is subtracted. The remaining 
signal is the flux generated by /, plus a residual, ellip- 
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FIG. 1: (Color online) <& a -I curves for a ring with w = 120 nm 
and R — 1.2 /im, fitted to a one-order-parameter Ginzburg- 
Landau model. The selected curves represent the regimes 
discussed in the text: (a) no transitions, (b),(c) hysteretic 
and (d) thermal equilibrium. 



tic sensor background which is negligible at lower T and 
unambiguously distinguishable from the ring response at 
higher T, where fluxoid transitions occur. Details on the 
technique will be given elsewhere 17 1. 

Rings with w < 120 nm show no two-OP effects. At 
low T [Fig. HJa)], there are no fluxoid transitions at the 
experimental time scale and field sweep amplitude. At 
higher T [Fig. [T^b), ( C )L we observe hysteretic transi- 
tions, which become non-hysteretic near T c [Fig. Old)]. 

We obtain the theoretical 4> a -J curve of an individual 
fluxoid state, n, directly from ID, 1-OP GL [H |: 
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where <p = $ a /$o, $o = h/2e, and d is the total film 
thickness, n is the phase winding number of the GL-OP 
i>(x) — \ip\e mx / B \ where x is the position along the ring's 
circumference, tp — n < i?/^GL gives the linear response 
of the London limit, while the cubic term arises from 
pair breaking. Because wd <C A 2 , the self inductance 
can be neglected. Although the width of some rings is 
several £gl at low T, the ID approximation is justified 
here because w <C R and H a <C H c2 [H| . 

Close to T c [Fig. [Tfd)], transitions are fast enough to 
model the experimental <fr a -I curves as a thermal average 
over all possible states, 
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with E n ((p) = -$ In dif'I n ((p'). We set £ GL = when 
substituting Eq. {l} into Eq. (|2]) because the thermal 
rounding dominates the cubic term [23| . The free param- 
eters in the fit are A -2 , three background parameters, a 
small offset in tp, the pickup loop-ring inductance M coup , 
and the field coil-ring inductance. The fitted inductances 



are consistent with less accurate geometrical estimates, 
and are used at lower T. 

The hysteretic curves [Fig. [ljb), (c)] are rounded 
from averaging over a distribution of thermally activated 
transitions [19( near a typical ip = <p t + n. We model 
them by combining Eq. ([T]) with occupation probabili- 
ties p n ((p — n) obtained from integrating the rate equation 
dp n /dt = -p n /T exp(-E act ((p - n)/k B T). For fitting, 
To and the activation energy E act are absorbed into 4>t- 
The free parameters are <pt, dE act /d(p(4>t), (,gl, A" 2 , and 
three background parameters. (p t increases very weakly 
with increasing field sweep frequency, as expected for 
thermally activated behavior. 

For T -c T c , where no transitions are seen, Eq. |T]) fits 
the data with three parameters corresponding to £gl, 
A -2 , and a constant background. 

The above models result in excellent fits for measured 
&a-I curves except at w — 190 nm as discussed below. 
The fitted values for A -2 and £gl are shown in Fig. [2] 
for rings typical of each w. All 14 measured rings with 
w < 120 nm show T dependences similar to the 1-OP 
phenomenological expressions A(T)~ 2 = A(0)~ 2 (1 — t 4 ) 
and £gl {T) = £ G L(0)V(l + * 2 )/(l-n with * = T/T c . 
For w > 190 nm, A _2 (T) has a high-temperature tail to 
an enhanced T c and a peak in ^gl near but below T c . 
Both effects are most pronounced at w = 190 nm. The 
$a-I curves remain hysteretic well into the tails, showing 
that the enhanced T c is not a fluctuation effect. 

The most striking features of the <& a -I curves for a 
W = 190 nm ring [Fig. [3fa)-(f)], typical for all six mea- 
sured rings with w = 190 nm and R > 0.8 ^m [2^, are 
transition points in the lower T hysteretic region that 
are not a function of if — n only, and reentrant hysteresis. 
The latter is qualitatively related to the local maximum 
in £ G l{T) [Fig. Gib)], since in 1-OP GL, a state be- 
comes unstable at \<p - n\ > y/R 2 /£ GL + 1/2/V3 
However, rather than an increase of the transitions point 
(pt upon raising T, as expected for a decreasing £,gl{T), 
the amplitude of those transitions is reduced until a non- 
hysteretic curve with a flattening too pronounced to be 
described by Eq. ([TJ) [Fig. [3Jd)] appears. With a further 
increase of T [Fig. [3je), (f)], this flattening also disap- 
pears and the $ a -7 curves evolve similarly to Fig. [T] This 
strongly suggests the existence of two OPs, one causing 
the small ripples in Fig. [31(c) , and one with a larger T c 
adding the large background response and causing the 
tail in A -2 (T). The additional phase winding number 
from a second, coupled OP also explains the irregular 
transitions in [Fig. (Ha)-(c)]. Although also breaking 
strict flux periodicity, finite line width corrections to ID, 
1-OP GL [18( would be much smaller and more regular. 
A vortex pinned in the annulus could in principle lead to 
similar $ a -/ curves, however w is too small compared to 
£gl to accommodate its core, and it cannot explain the 
T dependence of A -2 and 

A two-OP GL-model consisting of two standard ID 
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FIG. 2: (Color online) (a) A -2 and (b) £gl for rings repre- 
sentative of each annulus width w. The discrete symbols are 
obtained from fits to & a -I curves. Continuous curves repre- 
sent fits to phenomenological expressions (w = 120 nm) or 
the two-order-parameter GL model (w > 190 nm), which was 
only fitted above 1.2 K where GL applies. 



GL free energy functionals and a coupling term indeed 
reproduces the peculiar features of Fig. [2] and [3] 
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We define i/'i to have the lower T c . If both components 
have the same n, one can make the usual ansatz ipi(x) = 
\tj)j\e mx / R . Minimizing ^ with respect to and \^ 2 \ 
results in excellent fits to $ a -I curves as in Fig. [2Id). 
At small if, both OPs contribute significantly, but as ip 
is increased, pair breaking strongly reduces \ipi\, whose 
£gl diverges near T c> i. In the absence of the more stable 
ip2, "01 would undergo a fluxoid transition much before 
l^i | could be suppressed that much. 

We extracted effective values of A~ 2 (T) and £gl(T) 
from such modeled $ a -7 curves (and also from fits to 
datasets similar to Fig. Eld)) using a small tp expan- 
sion analogous to Eq. {T]). Assuming a linear T depen- 
dence of ai and a2, this procedure can reproduce the 
observed form of \~~ 2 (T) and £,gl{T) as demonstrated 
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FIG. 3: (Color online) (a)-(f): & a -I curves for a ring with w 
— 190 nm and R = 2 ^m, where the coupling between the two 
order parameters is weak. In (a) - (c), states with two different 
fluxoid numbers lead to the observed transition pattern. In 
(d), no transition occurs because one component is stabilized 
by the other. Dotted lines show the contributions of tpi and 
ij)2 derived from the model, (e) and (f ) , taken above the lower 
T c , reflect the response of a single order parameter, (g)-(i): 
Results from the two-OP model corresponding to (a)-(c). (j): 
Schematic of film structure (cross section). 



by the fits in Fig. [51 The knee in A _2 (T) corresponds 
to the lower T C] i, above which the amplitude of be- 
comes very small. From w = 190 nm to w = 370 nm, 
the coupling strength 7 increases by a factor 30-50, by 
far the largest line width dependence of all fit parameters 
[25( 1 . The resulting stronger proximitization smears out 
the two-OP features. It appears that ipi not only has 
a lower T c , but also a smaller l e than ip2- Fit parame- 
ters obtained from $ a -I curves and the curves in Fig. [5] 
are consistent within about a factor 2. This discrepancy 
may be due to fluctuations of ipi at large ip, or a nontriv- 
ial phase difference which are not considered in our 
model for Q a -I curves. 

To model the hysteretic <& a -I curves [Fig. [3](a)-(c)] , we 
calculate the activation energies E ac t(<p ~ n ) f° r transi- 
tions in either ipi or ijj 2 by numerically computing saddle 
point solutions of the GL - equations obtained from vari- 
ation of ©, assuming that the relevant saddle points 



4 



evolve continuously from those of the uncoupled system 
upon increasing 7. We then derive the complete $ a -I 
curve assuming that a transition in occurs whenever 
Eactiv — n)< KikgT, where the Ki 2 are treated as phe- 
nomenological parameters of order unity. Of several sim- 
ulation runs with various parameters similar to those ob- 
tained from the fits, the one shown in Fig. [3jg)-(i) gave 
the best similarity with the data [Fig. B^a)-(c)]. While 
the uncertainty of the fit parameters |25j and simplicity 
of the model forbid a more quantitative comparison, the 
simulations show that mctastable states with n% =/= n 2 
arc key to understanding the observed § a -I curves. Due 
to the coupling between the OPs, the variation of the 
relative phase is soliton-like 0], similar to a Josephson 
vortex, however with the phase gradient along the junc- 
tion being mostly due to the kinetic rather than the mag- 
netic inductance. This soliton corresponds to the domain 
walls originating from unquantized vortices in bulk sam- 
ples Q. Intuitively, it is formed upon increasing ip when 
the larger £gl of ipi causes ipi to become unstable at a 
smaller ip than ip 2 > and the coupling is weak enough for 
-02 to stay in the same state. However, the soliton energy 
makes it less stable than states with rii = n 2 . This leads 
to the observed step-like transition sequence. 

Apart from the general agreement with the simulations 
and the irregular transitions, the most direct experimen- 
tal evidence for states with m 7^ n 2 arc branches of <f> a -I 
curves that are shifted relative to each other by less than 
one $0 horizontally, such as around <& a /$o ~ ±3 in Fig. 
[3jb). Individual, unaveraged field sweeps show that this 
is not an effect of averaging over different transition path- 
ways, contrary to the features at $ a /$o ~ ±0.3. 

The emergence of two OPs can be explained by the 
temporary drop of the deposition rate during the metal- 
ization. At the lower rate, more oxygen was co-deposited 
to form a tunneling barrier [see Fig. [Uj)]. Different 
oxygen concentrations and/or grain sizes led to different 
values of T c and £gl in the two superconducting layers. 
It is known that PMMA outgases significantly and that 
thinner lines are affected more [2l[. This likely caused 
the line width dependence of the coupling and the com- 
plete oxidization of one of the superconducting layers for 
w < 120 nm, where we found no evidence for two OPs. 
The critical Josephson current densities estimated from 
the inferred values of 7 support this picture [25j ] . 

This fabrication result was unintentional but fortu- 
itous. While the parameters should be tunable a pri- 
ori with controlled exposure to oxygen gas, it would be 
difficult to obtain reproducible results from outgasing re- 
sist. Nevertheless, the data and analysis draw a clear 
picture of a two-OP superconductor with GL parameters 
that depend consistently on w. This dependence and the 
occurrence of two different T c 's allowed the study of a 
wide range of parameters. The insight thus gained may 
be used to design similar experiments on intrinsic two- 
component superconductors. The creation and detection 



of h/4e vortices in Sr2Ru04 would be of particular inter- 
est 0, Q . Since their energy is logarithmically or linearly 
divergent in the sample size @,|3], they might only be 
accessible as metastable states in mesoscopic samples, 
similar to the soliton states discussed here. 

In conclusion, we have explored effects emerging from 
two coupled order parameters in mesoscopic supercon- 
ducting rings. The most interesting ones are an anoma- 
lous temperature dependence of the average superfluid 
density A -2 and effective coherence length £gl, and a 
qualitative modification of the behavior of phase slips 
related to previously predicted metastable states with 
two different phase winding numbers 0, 0| and a soliton- 
shaped phase difference 

This work was supported by NSF Grants No. 
DMR-0507931, DMR-0216470, ECS-0210877 and PHY- 
0425897 and by the Packard Foundation. Work was per- 
formed in part at the Stanford Nanofabrication Facility, 
which is supported by NSF Grant No. ECS-9731293, its 
lab members, and industrial affiliates. We would like to 
thank Per Delsing, Egor Babaev and Mac Beasley for 
useful discussions. 



* Electronic address: hcndrikb@stanford.edu 
[1] C. Buzea and T. Yamashita, Superconductor Science & 

Technology 14, R115 (2001). 
[2] E. Babaev, Phys. Rev. Lett. 89, 067001 (2002). 
[3] E. Babaev, Nuclear Physics B B686, 397 (2004). 
[4] Y. Tanaka, Phys. Rev. Lett. 88, 017002 (2002), note that 

there appears to be an error in the boundary conditions, 

which affects some of the conclusions. 
[5] A. Gurevich and V. M. Vinokur, Phys. Rev. Lett. 97, 

137003 (2006). 
[6] A. Gurevich, Phys. Rev. B 67, 184515 (2003). 
[7] C. Noce and L. Maritato, Phys. Rev. B 40, 734 (1989). 
[8] M. Stone and S.-B. Chung, Phys. Rev. B 73, 014505 

(2006). 

[9] S. Das Sarma, C. Nayak, and S. Tewari, Phys. Rev. B 
73, 220502(R) (2006). 
[10] E. Babaev, A. Sudbo, and N. W. Ashcroft, Nature 431, 
666 (2004). 

[11] P. G. Bjornsson, B. W. Gardner, J. R. Kirtley, and K. A. 

Moler, Rev. Sci. Inst. 72, 4153 (2001). 
[12] P. G. D. Gennes, Rev. Mod. Phys. 36, 225 (1964). 
[13] R. Meservey and B. B. Schwartz, Equilibrium properties: 

comparison of experimental results with predictions of the 

BCS theory. (Marcel Dekker Inc., New York, 1969), pp. 

117 - 84. 

[14] R. B. Pettit and J. Silcox, Phys. Rev. B 13, 2865 (1976). 
[15] R. W. Cohen and B. Abeles, Phys. Rev. 168, 444 (1968). 
[16] M. E. Huber et al., (unpublished). 

[17] N. C. Koshnick, H. Bluhm, M. E. Huber, and K. A. 

Moler, (unpublished). 
[18] X. Zhang and J. C. Price, Phys. Rev. B 55, 3128 (1997). 
[19] J. S. Langer and V. Ambegaokar, Phys. Rev. 164, 498 

(1967). 

[20] D. Y. Vodolazov, F. M. Peeters, S. V. Dubonos, and A. K. 



5 



Geim, Phy. Rev. B 67, 054506 (2003). 
[21] P. Dubos, P. Charlat, T. Crozes, P. Paniez, and B. Pan- 

netier, J. Vac. Sci. Technol. B 18, 122 (2000). 
[22] F. von Oppen and E. K. Riedel, Phys. Rev. B 46, 3203 

(1992). 

[23] Fluctuation effects within a few mK of T c [22l |. where it 
would be inadequate to set £ql = 0, will be discussed in 



Ref. [I3|. 

[24] For rings with R = 0.5 fim, 4> a was too small to observe 
those features. 

[25] See auxiliary document at 

http : //www . st anf ord. edu/group/moler /publications . html 
for a discussion of the fit procedures and parameters. 



